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The continuum variability of MCG-6-30-15: A detailed analysis of 
the long 1999 ASCA observation 
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ABSTRACT 

We report on an analysis in the 3-10 keV X-ray band of the long 1999 ASCA observation 
of MCG-6-30-15. The time-averaged broad iron K line is well-described by disk emission 
near a Schwarzschild black hole, confirming the results of earlier analyses on the ASCA 1994 
and 1997 data. The time-resolved iron-line profile is remarkably stable over a factor of three 
change in source flux, and the line and continuum fluxes are uncorrelated. Detailed fits to the 
variable iron-line profile suggest that the active region (parametrized by the best-fit inner and 
outer radii of the accretion disk) responsible for iron line emission actually narrows with in- 
creasing flux to a region around 4-5 r g . In contrast to the iron line, the power-law continuum 
exhibits significant variability during the 1999 observation. Time-resolved spectral analysis 
reveals a new feature in the well-known photon index (T) vs. flux correlation: V appears to 
approach a limiting value of T ~ 2.1 at high flux. Two models are proposed to explain both 
the new feature in the V vs. flux correlation and the uncorrelated iron-line flux: a phenomeno- 
logical two power-law model, and the recently proposed "thundercloud" model of Merloni & 
Fabian (2001). Both models are capable of reproducing the data well, but because they are 
poorly constrained by the observed T vs. flux relation, they cannot at present be tested mean- 
ingfully by the data. The various implications and the physical interpretation of these models 
are discussed. 
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1 INTRODUCTION 

In recent years, a self-consistent and robust picture has gradually 
emerged in order to explain the most important features of AGN 
X-ray spectra. In the standard model, a central supermassive black 
hole accretes matter in the form of a cold, optically-thick disk. A 
hot, optically-thin corona above and below the accretion disk in- 
verse Compton scatters soft optical/UV photons from the disk, pro- 
duci ng a power-law continuum in the X-ray band (Z dziarski et al. 
1994). At the same time, the corona illuminates the disk, giving 
rise to a reflection hum p at higher energies and a prominent iron K 
line at around 6.4 keV (|Guilbert & Rees 1988j; L ightman & White 
1988). 

Extensive observations of MCG-6-30-15, a bright, nearby 
(z — 0.0078) Seyfert 1, have been especially useful in the study 
of AGN X-ray spectra. In fact, MCG-6-30-15 is perhaps the most 
well-studied AGN in the X-ray band, owing to its brightness, strong 
iron-line, and extreme variability. Because both time-averaged and 
time-resolved spectral analysis are feasible for this source, its phys- 
ical properties can be probed in great detail. The many interesting 
features of the time-averaged spectrum, such as the broad iron line, 
the reflection continuum, and the warm absorber, have proven to be 
extremely useful in constraining the geometry and physics of AGN. 

The study of spectral variability can provide another win- 



dow into the interesting properties of AGN, although it is gen- 
erally much more difficult and ambiguous than the simpler task 
of describing the time-averaged spectrum. Correlations (and lack 
thereof) have been observed between many different time-resolved 
spectral parameters, and there is much debate over the theoreti- 
cal interpretation. To mention two relevant examples, analyses by 
Lee et al. (2000) and Vaughan & Edelson (2001) of a long 1997 
RXTE dataset detected a steepening of the continuum spectrum 
with increasing flux in MCG-6-30-15. Meanwhile the iron line 
flux, though variable, did not appear to be correlated with the con- 
tinuum flux. 

In this paper, we present an analysis of the ASCA 1999 long 
observation of MCG-6-30-15 in the 3-10 keV X-ray band, inves- 
tigating both the time-averaged properties of the spectrum and the 
time-resolved spectral variability. Section |2| briefly describes the 
details of the 1999 observation and the status of the detector at the 
time. In Section^ the analysis of the time-averaged spectrum is out- 
lined. The time-averaged spectrum shows the presence of a strong, 
relativistically-broadened iron line, confirming the results of earlier 
analyses on the 1994 and 1997 ASCA data. Section ^ contains the 
results of the variability studies. The T vs. flux relation is inves- 
tigated, and a break in the usual correlation is observed. Iron-line 
variability with respect to flux is also studied, and in particular, the 
iron-line and continuum fluxes are found to be uncorrelated. Possi- 
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ble interpretations of the flattening of V are considered in Section 
including an application of the recently proposed "thundercloud 
model" of Merloni & Fabian (2001). Finally, Section ^ summarizes 
the results of the analysis. 



2 OBSERVATIONS AND DATA REDUCTION 

MCG-6-30-15 was observed with ASCA from 1999 July 19 to 
1999 July 29. The total integration time was ~ 910 ks. The data 
were filtered and spectral files were prepared following the method 
of Iwasawa et al. (1999). The Solid-State Imaging Spectrometer 
(SISO and SIS1) was operated in Faint mode for the duration of the 
observation, and the Gas Imaging Spectrometer (GIS2 and GIS3) 
was operated in PH mode. Good exposure time amounted to ap- 
proximately 440 ks for the SIS and 360 ks for the GIS. The FWHM 
energy resolutions at 6.4 keV for the SIS and GIS were ~ 350 eV 
and ~ 490 eV, respectively. This represents a degradation of a fac- 
tor of 3-<\ in SIS energy resolution from launch, resulting from ra- 
diation damage to the CCD chips. The GIS detectors, on the other 
hand, being gas-based, are relatively unaffected by the effects of 
high energy radiation, and so their spectral resolution has remained 
very stable over time. 
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Figure 1. The broad iron Ka line from the 1999 long ASCA observa- 
tion of MCG-6-30-15. The line has been constructed from a power-law 
plus reflection (pexrav) fit to the underlying continuum. Plotted here is 
the total time-averaged SIS data, corrected for detector efficiency. Note 
the extremely broad, two-horned profile - the FWHM of the blue core 
alone is ~ 2-3 times the detector resolution. The solid line shows the 
best-fit diskline model for a Schwarzschild black hole, with x 2 = 
1885.8/1889 dof for the complete pexrav plus diskline model fit 
to SIS+GIS. The best-fit diskline model parameters are given in Table 

s 



3 TIME-AVERAGED SPECTRAL FEATURES 

The total integrated count rates between 3-10 keV for the current 
observation were 0.540 and 0.518 cts/s for the SIS and GIS, re- 
spectively. The corresponding time-averaged flux in the 3-10 keV 
band was 3.20 x 10 -11 erg cm -2 s _1 . This represents a 20% in- 
crease over the 3-10 keV flux in the previous 1997 observation, 
indicating a significant brightening of the source between the two 
observations. 

The 3-10 keV spectrum of MCG-6-30-15 was described us- 
ing two XSPEC model components: pexrav for the power-law 
continuum and coldj-eflertkmfrorrHhe_ara 
& Zdziarski 1995); and diskline for the relativistically- 
broadened iron Ka line at a rest-frame energy of 6.4 keV from 
an accr etion disk around a Schwarzschild black hole (F abian et 
al. 1989). In addition, pexr av was relativistically blurred using 
rdblur (Fabian et al. 1989), and the entire model was modified 



by cold Galactic absorption of 7V H _ 4.06 x 10 cm (E lvis et 
al. 1989). 

The free parameters consisted of the inner and outer disk ra- 
dius, disk inclination angle (diskline and rdblur); the total 
line flux (diskline); and the power-law slope and normaliza- 
tion (defined to be the photon flux at 1 keV of the primary power 
law) (pexrav). The inclination angle relative to the observer of 
the reflecting material in pexrav was tied to the corresponding 
diskline parameter. Since the other pexrav parameters affect 
primarily the high energy spectrum, which is inaccessible to ASCA, 
they could not be constrained in this analysis, and so they were 
kept fixed at value s consistent with recent measuremen ts by Bep- 
poSAX and RXTE ( fcuainazzi et al. 1999j |Lee et al. 199$ ). Thus the 
pexrav cut-off energy was fixed at 150 keV, the reflection frac- 
tion was fixed at Q/2tt = 1 (where Q. is the solid angle subtended 
by the reflecting material as seen by the primary source), and the 
iron abundance was fixed at 1 Zq. 

The model was fit simultaneously to the data from all four 
detectors between 3-10 keV. A good fit was obtained with \ 2 — 



1885.8 for 1889 degrees of freedom. The best-fit continuum power- 
law slope and normalization are T — 2.082 ± 0.009 and A^-io = 
(1.90±0.03) x 10~ 2 ph s _1 cm~ 2 . The power-law slope is slightly 
higher than the s lope r 97 = 1.94+° n® measured for the 1997 (I wa- 
sawa et al. 1999) dataset. The difference in power-law slopes re- 
mains at the 2-a level even after systematic differences between the 
two analyses, such as the modification of pexrav with rdblur 
in this analysis and slightly different fitting procedures, are taken 
into account. 

The best-fit parameters for the diskline model are shown 
in the first row of Table p], and the efficiency-corrected iron line 
profile is shown in Fig. [l]. Generally speaking, the diskline pa- 
rameters agree quite well with the parameters measured in the 1997 
observation, confirming not only the presence of a broad iron K line 
in MCG-6-30-15, but also the properties of the accretion disk re- 
sponsible for the line emission. It is interesting to note, however, 
that while the 1997 data preferred an inner radius of 6.7r g , here 
ri n pegs at 6r g , the radius of minimum stability for a non-rotating 
black hole and the lower bound in the diskline model. This, 
and the relatively high value of a, suggest the presence of emission 
from within 6r g , which is generally thought to be the signature of 
a rotating black hole. 

To check for the significance of any emission from within 6r g , 
a model for a maximally spinning (a — 0.998) black hole was also 
fit to the data between 3-10 keV. The model is completely analo- 
gous to the one described above for a Schwarzschild black hole, 
with diskline and rdblur replaced with the XSPEC model 
components laor and kdblur, which describe line emission and 
relativistic blurring from an accretion disk around an extreme Kerr 
black hole (Kojima 1991; ^aor 1991). Again, we obtained an ex- 
cellent fit with T = 2.096 ± 0.009 and \ 2 = 1858/1889 dof. The 
best-fit laor parameters are shown in the second row of Table 
|IJ As indicated by the equally acceptable values of x 2 , the time- 
averaged line profile cannot constrain the spin of MCG-6-30-15 
given the current quality of the data. It is worthwhile to note, how- 
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Table 1. The best-fit parameters of the diskline and laor models for the iron-line profile of the 1999 ASCA observation of MCG-6-30-15. (1) power-law 
index of the disk radial emissivity profile (e oc r~ a ); (2) inner radius of the disk; (3) outer radius of the disk; (4) disk inclination angle; (5) iron-line intensity; 
(6) iron-line equivalent width. All quoted error bars are at the 1-cr level and include variations in the underlying continuum. 



ever, that in the best-fit Laor model, the entire line-emitting region 
of the accretion disk shifts inward and stretches all the way down 
to ~ 2 gravitational radii, resulting in a smaller emissivity index, 
since it is no longer necessary to concentrate as much emission so 
close to the black hole. 



4 TIME-RESOLVED SPECTRAL FEATURES: 
VARIABILITY STUDIES 

The background- subtracted SIS 3-10 keV light curve for the 1999 
observation is shown in Fig. |^ with 256-s bins. As in 1997, the 
source appears to be highly variable, with the 3-10 keV count rate 
changing in the brightest flares by factors of 2-3 on timescales of 
less than a day. 

4.1 Hardness Ratios 

Hardness ratios provide a straightforward, "zeroth-order" means of 
investigating spectral variability. Being model-independent, they 
can provide direct insight into overall time-resolved spectral fea- 
tures. Also, because they do not require any spectral fitting, hard- 
ness ratios are free of any biases introduced by a specific fitting 
procedure. A total of three hardness ratios were examined: 
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The hardness ratios were calculated in each orbital period (~ 5700 
s) using the total SIS+GIS count rate in the relevant energy bands. 
They are shown in Fig. |3| plotted against count rate. Ri and R2 
show a definite negative correlation with count rate (x 2 = 240.7 
and 242.8, respectively, for 158 dof against a constant hypothe- 
sis), indicating that the spectrum steepens as the source becomes 
brighter. On the other hand, Rs does not appear to show much ev- 
idence for this correlation, although we note that any such corre- 
lation could be obscured by scatter due to the significantly poorer 
signal-to-noise of the Rz band. 

The hardness ratio trends shown here agree well with those 
found by Lee et al. (2000) in their recent analysis of the 1997 RXTE 
observation, suggesting that the observed correlation between spec- 
tral steepness and source brightness is an intrinsic and general fea- 
ture of MCG-6-30-15. Since we expect the spectrum between 3- 
10 keV to be well-modeled by a simple power-law plus a relativisti- 
cally broadened iron-line, the observed hardness ratio trends should 



correspond directly to a positive continuum index vs. count rate 
correlation. As we shall see in the next section, this is almost what 
is observed. 



4.2 The 12 Count-Rate Bins 

The study of correlations between various spectral parameters and 
source flux was carried out using twelve bins defined using hor- 
izontal slices of the 3-10 keV light curve of the GIS3 detector. 
The procedure was as follows: the 3-10 keV GIS3 light curve was 
binned into orbital periods (6 x 10 3 s). From this light curve, twelve 
horizontal slices in count rate were chosen such that each slice con- 
tained approximately the same number of counts. These twelve hor- 
izontal slices are shown in Fig. ^ along with the GIS3 light curve, 
binned into orbital periods. The horizontal slices were then con- 
verted into timing filters (i.e. a union of disconnected vertical time 
slices, each with GIS3 count rate within the specified range). Fi- 
nally, the timing filters were applied to the data from the other three 
detectors (SIS0, SIS 1 , GIS2) to produce the desired spectral files. 

We decided to use horizontally sliced count-rate bins instead 
of the more customary vertically sliced time bins for several rea- 
sons. First, binning in count rate allowed us to increase signal-to- 
noise while preserving the large amount of variability observed in 
the source. Time bins of comparable size would have smeared out 
most of the interesting flares and minima, thereby destroying the 
greater part of the correlations to be described below. Secondly, us- 
ing count-rate bins in some sense averages out the dependence of T 
and other spectral parameters on properties of the source other than 
count rate and focuses on the sole determination of these spectral 
parameters as a function of count rate. 



4.3 Photon Index vs. Count Rate Correlation 

The first correlation we investigated was the standard relation be- 
tween the continuum photon index (r) and count rate. A positive 
correlation between F and source flux is believed to be a feature 
common to most Seyfert l's (Perola et al. 1986, Nandra et al. 1991, 
Ptak et al. 1994; see also the review by Mushotzky et al. 1993). 
Recent studies using high-quality data have firmly established the 
presen ce of a correlation in a num ber of AGN, including MCG-6 - 
30-15 dVaughan & Edelson 200 1| ), NGC 5548 ( JChiang et al. 2000| ), 
NGC 7496 ( |Nandra et al. 2000| ), and IC 4329A ( pone et al. 2000 ). 
The 1999 ASCA observation of MCG-6-30-15 is the longest and 
highest signal-to-noise observation of this source, and as such it of- 
fers an unprecedented opportunity to determine the photon index as 
a function of flux. However, the result of our analysis is surprising. 
Instead of finding a simple positive correlation between T and flux 
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Figure 2. The ASCA SIS light curve of MCG-6-30-15 in the 3-10 keV band. The data has been background-subtracted and is shown here in 256-s bins. The 
source was observed by ASCA from 19 July 1999 to 29 July 1999 for a total of 910 ks. 



as expected from previous analyses, we report here a new trend: a 
flattening of Y at high count rate. 

r was determined using the blurred power-law plus reflection 
model described in Section H. For simplicity, the iron line was not 
simultaneously determined in these fits, so as to avoid any model 
degeneracies and biases (e.g. between the power-law slope and the 
iron-line flux) that might arise due to the reduced signal-to-noise 
of the time-resolved data. As in the time-averaged fit, the free pa- 
rameters were T and the continuum normalization. Fig. ^| shows 
the best-fit photon indices plotted vs. the average 3-10 keV count 
rate of each flux bin. Absorption (warm or otherwise) is not ex- 
pected to play a significant role above 3 keV in this source (al- 
though we test this assumption below), and X-ray reprocessing fea- 
tures such as the iron line and the reflection hump are estimated to 
contribute at most 15% to the total 3-10 keV flux. Therefore the 
3-10 keV count rate should accurately reflect changes in the flux 
of the primary continuum. These changes are seen to result in a 
positive correlation of photon index with flux when the flux is low, 
in qualitative agreement with previous analyses of MCG-6-30-15 
and other sources. However at higher flux, the photon index grad- 
ually levels off to around T ~ 2.1-2.2. A x 2 -test against a linear 
correlation gives \ 2 = 15.3/10 degrees of freedom, correspond- 
ing to a rejection of the linear correlation at 88% confidence. The 
fit is dramatically improved with the power-law functional form 
r = T - KF~ S , which yields X 2 = 6.65/9 dof with a limiting 
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and a (negative) power-law index of 



We believe this to be the first indication in an AGN 



of a saturation of the photon index at high flux. Should this result 
prove to be real, it could have serious implications for theoretical 
models of coronal activity and continuum production, as we shall 
discuss in Section ^|. 

We have tested the saturation of T(F) at high flux and the 
result appears to be robust. Varying the fit prescription and frozen 
model parameters had little effect on the general trend. Similarly, 
it made little difference to the shape of the correlation whether rel- 



ativistic blurring of the continuum (rdblur) was included in the 
fit. 

Obvious from inspection of Fig. ^ is the important role played 
by the highest flux bin (J^-io ~ 1.6 cts/s) in determining the 
nature of the relation between T and flux. Excluding the high- 
est flux bin does improve the linear fit to a marginally acceptable 
X 2 = 11.5/9 dof. However, the residuals exhibit additional struc- 
ture that still suggest a saturating T. This is confirmed by the fact 
that fitting with the power-law functional form results in a best-fit 
nearly identical to the full 12-bin fit. The stability of the power-law 
fit, along with the fact we can think of no good physical reason why 
the highest flux bin should be excluded from the relation, give us 
confidence that the apparent saturation of V is not being caused by 
a single bin. 

Other effects, such as correlated errors and model-specific bi- 
ases, could be more subtle. MCG-6-30-15 is known to contain a 
highly variable warm absorber, with the absor ption dominated by 
O VII and O V III edges at 0.74 and 0.87 keV ( Fabian et al. 1994 ; 
Lee et al. 2001). While the warm absorber is not expected to play 
a significant role above 3 keV, we have nonetheless tested for the 
possibility that variable warm absorption could have an effect on 
the r vs. flux relation. Severe detector degradation below 2 keV 
due to radiation damage prevented detailed and realistic warm ab- 
sorber fits from being done in the 1999 observation. Instead, we at- 
tempted to model the effect of warm absorption above 3 keV using 
an additional cold absorption model component (wabs, Morrison 
& McCammon (1983)). The fits were redone with the column den- 
sity of this extra wabs left as a free parameter. The fits actually 
preferred a rather high column density of ~ 10 22 cm~ 2 , which re- 
sulted in slightly softer photon indices. However, since the amount 
of absorption was not correlated with flux, the resulting V vs. flux 
relation was largely unchanged and still showed the saturation of T 
at high flux. 

A systematic study was conducted to ascertain whether there 



The continuum variability ofMCG-6-30-15 5 




0.6 

F +F 

3-4 4-7 



K 



o 

o 



00 




o 2 



2x10" 



4x10° 6xlO J 
Time (s) 



CR Bin 12 



3x10" 



Figure 4. Essentially the same light curve as in Fig. g, except derived from 
the GIS3 detector and binned with 6000-s bins. Shown here are the 12 hori- 
zontal slices of the GIS3 light curve that define the count-rate bins used for 
the spectral variability analysis described in this paper. 
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Figure 3. The hardness ratios R±, R2, and R3 of MCG-6-30-15, plotted in 
separate panels from top to bottom against the total count rate in the relevant 
energy bands. The clear downwards trend seen in the first two hardness 
ratios indicate a definite spectral steepening as the source becomes brighter. 



could be any count-rate dependent bias in the determination of 
10 simulated spectra were produced in each count-rate bin using 
the best-fit photon index and continuum normalization for that bin. 
Exposure times, response curves, and background spectra were de- 
rived from real data so as to mimic as closely as possible the actual 
fitting process. T was then measured for the 10 simulated spectra 

1 Strictly speaking, it is not the count rate, but rather the total number of 
counts (i.e. the signal-to-noise level) in each bin that determines the level 
of bias in the fit. Since the count-rate bins were defined so as to contain 
approximately the same number of counts, such a bias should in principle 
have a minimal effect on the shape of the trend. 



in each count-rate bin, and the average V was found to differ from 
the input photon index by less than 1% in every count-rate bin (well 
within the 1-cr errors on the actual T), with no systematic trend with 
respect to count rate. Thus the saturation of T at high flux cannot 
be accounted for by a flux-dependent bias in the fitting process. 

Finally, a similar analysis was carried out to measure the con- 
tinuum photon index as a function of flux during the 1997 ASCA 
observation of MCG-6-30-15. The result agrees remarkably well 
with that of the present observation, providing further evidence 
(along with the hardness ratios) that the T vs. flux correlation is a 
general feature of MCG-6-30-15. However, because the source was 
approximately 20% fainter in 1997, the bulk of F(F) for that obser- 
vation falls below the saturation regime of > 1 ct/s. Indeed, of the 
three long ASCA observations of MCG-6-30-15 in 1994, 1997, and 
1999, only the latter provides the exposure time and average source 
luminosity necessary to adequately probe the saturation regime of 
the r vs. flux relation. 



4.4 Iron-line Variability 

Temporal variations of the iron-line profile of MCG-6-30- 1 5 in the 
three long ASCA observations (1994, 1997, and 1999) were previ- 
ously studied by Matsumoto et al. (2001). They showed using the 
normalized RMS variability spectrum of the 1999 data that the vari- 
ability of the line band (5-6.6 keV) decreased faster than that of the 
continuum band as one moved to longer timescales. This suggests 
that the line varies on a shorter characterstic timescale than the con- 
tinuum. In addition, using double gaussian fits to the red and blue 
wings of the broad line, they found the line flux to be significantly 
variable during the 1994 and 1997 observations, but not during the 
1999 observation. The line equivalent width, however, was variable 
in all three. The equivalent widths of the red and blue wings showed 
a very weak positive correlation, with an unacceptably large scat- 
ter around a linear best-fit. The total equivalent width tended to be 
anti-correlated with continuum flux, contrary to the constant be- 
havior that would be expected from the simplest coronal disk-line 
models. 

In this analysis, we have studied in greater detail the relation- 
ship, if any, between properties of the iron line and the continuum 
flux. Detailed fits in the 3-10 keV band using the Laor disk-line 
plus continuum model described in Section ^| were performed for 
each of the twelve count-rate bins. The free parameters were the 
photon index, continuum normalization, inner and outer disk ra- 
dius, and total line flux. The inclination angle of the disk was frozen 
at its time-averaged best-fit value (see Table [j]), and the rest-energy 
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Figure 5. The continuum photon index as a function of the 3-10 keV count 
rate. All errors shown are l-cr. The photon index appears to saturate at high 
flux to r ~ 2.1, a new feature of the T vs. flux relation. 



of the line was fixed at 6.4 keV. The disk emissivity index was also 
frozen at its time-averaged value, since simultaneous determination 
of the emissivity index and the disk radii was subject to large de- 
generacies due to the reduced signal-to-noise of the time-resolved 
spectra. 

Fig. ^| shows the efficiency-corrected iron line profiles ob- 
tained in each of the 12 count-rate bins, along with the best-fitting 
model in each case. The energy and flux scales are the same in each 
panel. Even with the reduced signal-to-noise of the time-resolved 
spectra, a very broad iron line can still be clearly resolved in nearly 
every bin. In most cases, an extended red tail is seen stretching all 
the way down to (or even below) 3 keV, suggesting the presence of 
line-emitting matter extremely close (r ~ a few r g ) to the central 
black hole at almost all levels of source brightness. Overall, the line 
appears to be remarkably stable across the flux bins, with the ex- 
ception of perhaps the highest-flux bin, where the red wing seems 
to weaken considerably while the narrower blue wing remains at ~ 
6keV. 

The spectral parameters shown in Fig. [7] offer a more precise 
picture of the variability of the line with flux. The continuum pho- 
ton index F has been plotted in addition to the inner/outer radii 
and the line flux in order to demonstrate that the T vs. flux rela- 



tion found in Section 4.3 holds even when the iron-line and con- 
tinuum are simultaneously fit. Regarding the line parameters, we 
see that they accurately reflect the apparent stability that was seen 
in the line profiles of Fig. (excepting the highest-flux bin). At 
every flux level, n n is within 6r g , the radius of marginal stability 
for a Schwarzschild black hole. More importantly in the context of 
coronal models, the fits show that the iron-line flux clearly does not 
increase with the continuum flux. The iron-line flux is consistent 
with a constant level (x 2 = 7.2 for 11 d.o.f.) over nearly a fac- 
tor of three change in source count rate. The lack of any positive 
correlation between line and continuum flux appears to be robust, 
confirming the previous findings of Lee et al. (2000) and Vaughan 
& Edelson (2001) for MCG-6-30-15. 

The uncorrelated line and continuum fluxes present a serious 
problem for the simplest coronal models in which the iron line is 
produced by reflection of the primary continuum off a cold accre- 
tion disk, with the continuum generated in a corona directly above 
the disk. In such a scenario the iron line flux would be expected 
to respond proportionately to variations in the primary continuum, 
with the equivalent width remaining roughly constant with respect 



to the continuum flux. Our findings clearly indicate that the equiv- 
alent width of the iron line actually decreases with flux in MCG-6- 
30-15. Although a positive correlation between line flux and con- 
tinuum flux has been observed in some AGN, such as NGC 7314 



(Yaqoob et al. 1996), it appears to be completely lacking in other 



AGN besides MCG-6-30-15, most notably in NGC 3516 (N andra 
& George 1999) and NGC 5548 ( phiang et al. 200C| ). The fact that 
the line and continuum flux are not correlated in some AGN in- 
dicates that the reflection process, at least in these AGN, is more 
complex than that described by a simple coronal model. 

Stronger evidence for variation with flux is found for the outer 
radius r out . Since the disk inclination angle has been fixed to its 
best-fit time-averaged value, r out is relatively well-constrained by 
the well-defined peak and maximum energy of the blue wing. A 
constant level is a poor description of r out vs. count rate, with 
X 2 = 27.3/11 dof. The fit improves significantly when one fits 
with a linear relation (x 2 = 16.4/10 dof, best-fit slope — 4.4±1.3). 
Excluding the lowest-flux bin leaves the quality of the constant 
fit virtually unchanged, while the linear fit improves further to 
X 2 = 4.8/9 dof with slope -8.9 ± 1.9. Thus, overall, the outer 
radius appears to exhibit a significant anticorrelation with source 
flux, suggesting that the illumination pattern on the disk may shift 
closer to the black hole when the source brightens. It is intriguing 
that ri n also appears to exhibit a slight positive correlation with 
flux, although it is not statistically significant. The trend, if present, 
along with the anticorrelation of r out would imply that the line 
emission originates in a more localized part of the disk when the 
source is bright. One might expect this sort of behavior if variations 
in source flux were due to flaring active regions on the disk. The 
overall trends in r- m and r ou t could then be understood by assum- 
ing that the dominant active region was relatively stable throughout 
the 1999 observation and was located around 4-5 r g . 

A detailed, exhaustive analysis of the many possible degen- 
eracies between various fit parameters is beyond the scope of this 
paper. Thus we cannot rule out the possibility that the trends shown 
here are simply a result of correlated errors. Nevertheless, we were 
able to rule out one of the more obvious possible correlated errors 
with the following simple check. The plot of line flux vs. source 
flux in Fig. [7] strongly suggests an anti-correlation between line 
and source flux. Suppose however that the line flux was actually 
constant. Then since the count-rate bins contain approximately the 
same number of continuum counts, the higher count-rate bins, be- 
ing shorter in duration, would contain less counts from the iron 
line. Accurate determination of the line parameters would be more 
difficult, and in particular, the wings of the broadened line would 
be harder to detect. This would bias the line flux downwards and 
in addition spuriously narrow the line emitting region. To check 
that this systematic error due to decreased line counts is not the 
source of the behavior shown in Fig. ^| we have binned together 
the last two count-rate bins, and the spectral fits have been redone 
for this summed spectrum. The spectral parameters obtained for the 
summed spectrum follow the general trend, and they are found to 
lie between those of the two last count-rate bins, indicating that at 
least part of the variation in spectral parameters cannot be explained 
by a simple decrease in line counts, with the line flux constant. 



5 INTERPRETING THE T VS. FLUX CORRELATION 
5.1 The Two Power-Law Model 

We propose two different models that can simultaneously account 
for the observed saturation of the photon index at high flux and the 
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Figure 6. The broad iron line in twelve count rate bins, along with the best-fitting Laor model. Filled circles (open squares) correspond to data from the SIS 
(GIS). The bins are numbered in order of increasing flux. 



lack of correlation between the iron-line and continuum fluxes. The 
first is an ad hoc phenomenological model consisting of two power- 
law continua, one constant with flux and photon index (Fi,Ti), 
and the other variable with (^2^2). To test this model, a coarse 
grid in Ti and T2 was produced with Ti ranging from 1.75 to 
2.0 and Y2 from 2.2 to 2.3. F\ was fixed at 0.54 cts/s (chosen to 
correspond approximately to the minimum of the 3-10 keV light 
curve), because Fi and F\ proved to be highly degenerate param- 
eters. Simulated spectra for each value of Y\ and Y2 were pro- 
duced in every count-rate bin and were folded through the expo- 
sure times, response matrices, and effective areas of the actual de- 
tectors. Normalizations were determined in each bin by requiring 
the overall 3-10 keV count-rates of the two power laws to equal 
F\ and F2 = F — F\ respectively, where F was the average count 
rate of the actual data in that bin. The simulated spectra were fit 
with a single power-law model between 3-10 keV to produce V vs. 
flux relations for each value of Ti and T2. Shown in the left panel 
of Fig. HI along with the data is the "best-fit" for our coarse grid 
(Fi = 1.85, T 2 = 2.25 with F\ = 0.54; x 2 = 10.5/10 dof), 
interpolated using the functional form r = To — KF~ S . The term 
"best-fit" is slightly misleading, however: for every value of F2, 
there was a value of Ti that yielded an acceptable fit. In any event, 
the two power-law model, although poorly constrained by the data, 
is clearly capable of reproducing well the observed form of F(F). 

Aside from the fact that it models well the measured T vs. 
flux relation, the two power-law model is attractive because it of- 
fers a chance of decoupling the iron-line and continuum variability. 



Previous analyses, and the detailed results of Section %A , indicate 



strongly that the iron line flux is uncorrelated with the continuum 
flux in MCG-6-30-15. As we discussed above, this is very puz- 
zling in the context of simple coronal models in which the line is 
expected to respond directly to changes in the primary continuum. 
The two power-law model offers a potential solution to this prob- 
lem: by associating the iron line with reprocessed radiation from 
the constant component alone, the continuum is free to vary in re- 
sponse to the variable component without having any effect on the 
iron line. Althou gh the iron line is known to exhibit significant tem- 
poral variability (Iwasawa et al. 1999; Matsumoto et al. 2001), we 



see from Section 1.4 that it is essentially stable across flux bins. 
Therefore we must consider the possibility that the "constant" F\ 
component varies to some extent as well, but that its variations are 
smeared out in the ~ 80 ks count-rate bins. This would not change 
the r vs. flux relation predicted by the model. 

At present the two power-law model is basically a phe- 
nomenological model, with very little theoretical motivation. We 
suggest a few possible explanations for the origin of the two power- 
law components. One possibility is that the "variable" F2 power- 
law component originates in the inner edge of the accretion disk, 
where the solid angle subtended by the reflecting material is the 
least; while the "constant" F\ power-law is produced farther out 
(but still in the inner part of the disk). F2 could here be associated 
with magneti zed accretion jus t inside the marginally stable orbit 
( Krolik 1999 ), as simulations ( Hawley & Krolik 2001 ) show that 
this process can be very noisy. Another possible scenario is that 
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The best-fit parameters (continuum photon index, inner/outer disk radius, iron-line flux) of the iron-line variability fits, plotted vs. 3-10 keV count 



flares occurring at different heights above the corona are responsi- 
ble for the two power laws. DiMatteo et al. (1999) showed that the 
spectral properties of the soft-state and hard-state of the GBHC GX 
339-4 could be qualitatively explained by flares occuring close to 
and far above the accretion disk, respectively. Moreover, they found 
that the ionization state of the disk was quite different depending on 
the height of the flares, with flares closer to the disk naturally ioniz- 
ing it more. Therefore associating the F\ and F2 components with 
flares high above and close to the disk, respectively, would not only 
account for the difference in continuum shape, but also, if the flares 
close to the disk completely ionized away the iron-line, would ex- 
plain the lack of an iron line accompanying F2 . 

Finally, there is the possibility that a small but non-negligible 
non-thermal component is present in the corona. A hybrid 
thermal/non-thermal corona has been explored as a model for the 
soft-state Cygnus X-l X-ray spectrum and wa s found to explain 
well its overall shape (3ierlinski et al. 1999). Single scattering 



is that strong evidence for a thermal corona in MCG-6-30-15 is 
provi ded by observations of a cutoff in the hig h energy X-ray spec- 



of disk seed photons off the non-thermal component gives rise to 
a hard, high-energy power law, while multiple scattering off the 
thermal part produces a softer power-law tail above the thermal 
disk emission at lower energies. A similar situation in MCG-6- 
30-15 would lead to the harder F\ (softer F2) component of our 
two power-law model being associated with the non-thermal (ther- 
mal) part of the corona, and with the two states varying indepen- 
dently of one another. An obvious problem with this explanation 



trum (Guainazzi et al. 1999; Lee et al. 1999). We note, however, 



that even small non-thermal tails in the coronal electron distribu- 
tion can apparently have si zable effects on the shape and flux of 
the Comptonized spectrum ( Wardzinski and Zdziarski 2001 ). How 
well the observations limit the size of any non-thermal tails in the 
corona of MCG-6-30-15, as well as what affect these might have 
on the 3-10 keV X-ray spectrum, remains to be seen. 



5.2 The Thundercloud Model 

A more theoretically motivated model of coronal variability that 
predicts quite generally the saturation of the continuum photon in- 
dex at high flux was recently proposed by Merloni & Fabian (2001). 
Their so-called "thundercloud model" was based on the idea that 
active regions in the corona could, under the right conditions, trig- 
ger avalanches of neighboring flares, giving rise to progressively 
larger "magnetic thunderclouds" that are responsible for inverse 
Compton scattering soft disk photons into the observed X-ray con- 
tinuum. Larger active regions tend to be more luminous (since they 
contain more magnetic reconnection sites) and produce softer spec- 
tra. Therefore when the source is bright, the continuum is domi- 
nated by soft emission from large active regions. During the largest 
flares, the active regions approach the limiting case of a slab-like 
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Figure 8. Good fits to the data for the two power-law model (left) with Ti = 1.85, Y2 = 2.25, F\ = 0.54; and the thundercloud model (right) with r = 1.5, 
C = 0.001, and D = 1.4. For both models the predicted V vs. flux relations were interpolated using the function V{F) = T - KF~ S . As one can see 
from the plots here, it is impossible to rule out either model from the T vs. flux relation alone, although the thundercloud model certainly has the advantage of 
being more theoretically motivated. 



geometry covering the entire disk, resulting in a saturation of V at 
high flux. 

In the thundercloud model, the luminosity of an active region 
is assumed to scale with its size via the relation L(r) oc r D , where 
D is a free parameter of the model which describes both the internal 
structure of the region (a kind of fractal dimension) as well as any 
radial dependence there might be in the input radiation from the 
accretion disk. The number of active regions of size r generated at 
time t is chosen from a Poissonian distribution with mean n(r, t), 
which in turn is determined by requiring that the corona be in a 
stationary state. This gives the relation n(r) — n(r,t)to(r), where 
n(r) tx r~ p is total number of active regions of size r and to(r) is 
the typical lifetime of these regions. Here p is given by the relation 
p = 2D + 3 — 7, where 7 is the slope of the red noise part of the 
Power Density Spectrum (PDS). For comparison with MCG-6-30- 
15, the slope of the PDS was fixed at 7 _ 1.5, in agreement w ith 
the literature ( Nowak & Chiang 200C ; Vlatsumoto et al. 2001 ). In 
addition to D, the other two free parameters of the model were the 
covering fraction C, which determines the total number of active 
regions; and the optical depth r of the corona. The mass of the 
black hole, accretion rate, and maximum/minimum flare size were 
fixed at M BH = 10 7 M , m = 0.1, r min = 0.02, and r max = 4, 
respectively. These serve mainly to set the overall luminosity and 
time scales. 

The r vs. flux relation predicted by the thundercloud model 
was compared to the data using a grid with optical depth 0.5 < r < 
2, covering fraction 1 x 10 -4 < C < 2 x 10 -3 , and luminosity 
scaling index 0.6 < D < 1.6. In addition to the constraint given 
by the T vs. flux relation, the overall normalized RMS variance of 
the model light curve was required to lie within the range 20-30% 
to match the that of MCG-6-30-15. Shown in the right panel of 
Fig. ^ is a particularly good description of the data with r = 1.5, 
C = 0.001, and D = 1.4 (the RMS variance was 21%). Overall, 
the coronal optical depth was the most well-constrained parameter, 
being required to lie in the range 1 < r < 2. The constraint on r 
is a consequence of the fact that it determines the limiting value of 
the photon index: it is essentially the only remaining free parameter 
at the highest flux levels when the active regions limit to a disk- 
covering, slab-like geometry. 

It is worth considering briefly the physical implications of 
t > 1 for iron-line variability. Since a r > 1 corona tends to 



smother any reflected radiation from the disk, the line should es- 
sentially vanish during higher flux states, when the thunderclouds 
cover most of the accretion disk. Meanwhile, in lower flux states, 
the more or less continuous series of small, localized flares should 
lead to a roughly constant line. Referring back to the time-resolved 
line profiles of Fig. |^ and the iron-line vs. continuum flux plot of 
Fig. M, it is intriguing to see that, within the limitations of the poor 
signal-to-noise, the behavior of the iron-line profile as a function of 
flux does indeed appear to support this prediction of the thunder- 
cloud model. 

The covering fraction and scaling index were much more dif- 
ficult parameters to constrain, as they contribute principally to de- 
termining the overall curvature of T(F), which is poorly measured 
by the present dataset. We can, however, limit these parameters to 
C > 0.01% and D < 2. These bounds are related, since val- 
ues of C below 0.01% require D > 2 (and vice versa) in order 
to achieve realistic values of the RMS variance. This can be un- 
derstood by noting that when the total number of active regions is 
small, their luminosity must grow quickly with their size in order 
to produce significant variations in overall luminosity over time. 
Such large values of D in turn require p to be large, which gives 
a corona dominated by very small active regions. This results in a 
rather unrealistic light curve consisting of long periods of stability 
due to small active regions, punctuated with the occasional huge 
flare from a large active region. 



6 SUMMARY 

An analysis of the spectral variability of MCG-6-30-15 during the 
1999 long ASCA observation has been presented. The principal 
findings of the analysis will now be summarized. A model consist- 
ing of a power-law continuum, plus reflection and a relativistically 
broadened iron-line from the accretion disk, gave a good fit to the 
time-averaged 3-10 keV spectrum. The spectral parameters of the 
broad iron line were generally in good agreement with the results of 
previous analyses of MCG-6-30-15. While not statistically signifi- 
cant, the data seemed to prefer an inner radius of the disk within the 
radius of marginal stability r ms for a Schwarzschild black hole, in 
contrast to earlier observations which found a best fit inner radius 
just outside r ms . 
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A study of the time-resolved iron-line profile revealed a com- 
plete lack of any positive correlation between iron-line flux and 
continuum flux, confirming the results of previous analyses. In ad- 
dition, marginal evidence was found for a narrowing with flux of 
the line-emitting region on the accretion disk, suggesting that iron- 
line variability may be due in part to changes in the illumination 
pattern of the disk. 

The power-law continuum was found to exhibit a positive cor- 
relation between photon index and source flux that levelled off at 
high flux, with Y limiting around 2.1-2.2. Two different models 
were proposed to explain the observed saturation in Y, and both 
were capable of describing the trend well. The first model pro- 
posed was a simple phenomenological model consisting of two 
power laws, one constant with Y ~ 1.8, and the variable with 
r ~ 2.2. It was also suggested that the lack of correlation be- 
tween iron-line and continuum flux could be explained with this 
model by associating the line production with the harder, constant 
component. Possible origins for the two distinct power laws were 
considered, and coronal flares of varying heights above the disk 
and a small non-thermal component in the corona were put forth 
as viable hypotheses. The second model advanced to explain the Y 
vs. flux relation was the so-called "thundercloud model" of coro- 
nal variability recently proposed by Merloni & Fabian (2001). The 
thundercloud model also generally predicts a saturation in Y at high 
flux as the growing thunderclouds approach a slab-like geometry 
above the accretion disk. A detailed comparison with the thunder- 
cloud model constrained the optical depth of the corona to lie in the 
range r ~1— 2, and a lower limit on the coronal covering fraction 
of C > 0.01% was obtained. In addition, the thundercloud model 
with t ~1— 2 predicts a pattern of iron-line variability (roughly con- 
stant line profile when the continuum flux is low; decreasing line 
flux when the continuum flux is high) which accords well with the 
observed behavior of the line. 
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